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We model the emission spectrum of a quantum dot em- 
bedded in a (e.g. photonic crystal) nanocavity, using a ••••••••• 

semi-classical approach to describe the matter-field in- 

teraction. We start from the simple model of a quan- # # # # # # # # 

tum dot as a two-level system, and recover the result ex- m m ^ ^ ^ ^ ^ ^ 

pected from cavity quantum electrodynamics. Then, we w w w w 

study the influence of electron-acoustic-phonons interac- ••••••••• 

tion. We show that the surrounding semiconductor plays 

an essential role in the emission spectrum in strong cou- „ , ■ r , t , ^ , ■ . t,tt^ ■ , 

Schematic ot the LI cavity ot the 2D triangular PHC with em- 

P bedded quantum dot. 
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1 Introduction Semiconductor quantum dots (QDs), 
displaying a quantized spectrum of energy levels, can be 
considered as a solid state equivalent of atoms. Recently, 
they allowed to draw a parallel between semiconductor physics 
and Cavity Quantum Electrodynamics (CQED, formerly 
being restricted to the domain of atomic physics). This 
was made possible by tailoring the light-matter interac- 
tion in resonant photonic structures such as micropillars 
and photonic crystal (PHC) nanocavities. Indeed, various 
cavity related effects such as Purcell effect and Rabi split- 
ting have been experimentally observed ||T|. Understanding 
how the semiconductor environment surrounding a QD af- 
fects CQED properties is of great interest, as it gives a mea- 
sure of how the QD-cavity system departs from an ideal 
CQED system, thus assessing the feasibility of a QD im- 
plementation of quantum information technology. 

In this work, we model the emission properties of a 
QD embedded in a nanocavity (e.g. a LI defect in a PHC 
slab ifTOI ). We extend a semi classical formalism originally 
developed to study radiative properties of QDs in vacuum 
and planar resonators ||2| . The electromagnetic field is de- 
scribed in terms of Maxwell equations, and the dielectric 



structure of the nano-cavity is accounted for by means of 
the resolvent representation of the Green's function. A sus- 
ceptibility tensor models the QD degrees of freedom. In the 
simple two-level assumption, the well known CQED result 
is recovered. By including the electron-phonon coupling, 
we discuss how acoustic phonon sidebands 18J in the QD 
affect the strong coupling. 

2 Theory The classical electric field E {r, uj) obeys a 
usual Maxwell equation, where r is the position vector in 
3-D, e(r) is the spatially dependent dielectric constant and 
Xqd is the linear optical susceptibility tensor of the QD 
subsystem: 



V AV A£ {r,Lu) 



e{r)£ {r,Lu) 



-47r J dr'xQD [r, r' , u) ■ S (r', 



= 0. (1) 



While e(r) is a constant in the case of vacuum, it de- 
scribes in general the dielectric structure surrounding the 
QD. In the case of a PHC nanocavity - that we consider 
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here - it describes the spatial profile of the dielectric con- 
stant that varies between a low and a high index material. 

In order to obtain a Hermitian eigenvalue problem for 
the Maxwell equation, we define 



K. Then, for oj w we can write 



in) 



+ c{r,r',uj) . 



(9) 



and 



r = 




(2) 



(3) 



This leads to an equation for the hermitian linear dif- 
ferential operator T: 



r)Q(r,a;) 



(4) 



We are interested in modeling the linear response of the 
system to an input field So in the following way |l3|: 



£ir,u;) 



£o 



(5) 



with S (lu) the emission spectrum and ^ (r , w) a form fac- 
tor that gives the spatial and spectral dependence of the 
field emitted outside the cavity. One central quantity in 
our approach is the Green's function corresponding to the 
Maxwell equation (|4|i, defined as; 



Here, c(r , r') represents the contribution of all other modes, 
and is supposed to be small at a; « ojc- A complete nu- 
merical calculation of cavity eigenmodes ifTOl . would al- 
low to test this assumption. Alternatively, our approach can 
be justified in terms of the quasi-mode theory ifTTl . by as- 
suming weak coupling between an ideal undamped cav- 
ity mode and the vacuum electromagnetic field outside the 
cavity ifTTI . Then, the spectrum can be obtained from Eq. 
(|7]i and an appropriate description of the QD susceptibility 

iCQD (r, r', a;) = ^•f' (r, r) M^* (r', r') x (t^) , (10) 

where is the Bloch part of the interband dipole matrix 
element and !Z/(re,r;i) is the electron-hole-pair envelope 
wave function in the QD. 

3 CQED limit The CQED Umit describes a two level 
atom in a high finesse cavity Ul. This result can be trans- 
posed to the case of a solid state oscillator and a microcav- 
ity ||4|. We start from a single electron-heavy-hole transi- 
tion at cj = ti>o. with an additional non-radiative linewidth 
7o of the bare QD: 



g{r,r',Lu) = 5{r~r'). 



(6) 



1 



LJo - LU - 270 



(11) 



Inserting Eq. (|6j into Eq. ([Sj, it is straightforward to 
prove that the following integral equation holds: 

Q(r,w) = Qo ir,u;) 



(7) 



We assume that the cavity mode is much more extended 
than the electron-hole pair wave function, that are thus ap- 
proximated by dirac-deltas (tq is the position of the QD). 
We finally get 



-47r— / / dr'dr"g{r,r',u) 



Xqd {r',r",uj) 



where Qo is the solution for the bare cavity (i.e. the 
PHC nanocavity in the absence of an embedded QD). When 
the QD is embedded in the vacuum or in a planar microcav- 
ity, an analytical expression for the Green's function holds 
121 ■ In a PHC nanocavity this is in general not the case. Yet, 
a very effective simplification is provided by the use of the 
resolvent representation of the Green's function, in terms 
of a set of orthogonal eigenmodes <?„(r) of the cavity ||9l : 



i^<f 0(1-0) 



Atc(ujo - u) - 170) 



(12) 



(oJo - tAj - i^Qo){(^c - w - ifi) - g2 



with Aec = ej\/ — 1 - where ejv/ is the dielectric con- 
stant of the medium in which the cavity is formed, and Sb 
is an external field entering the system. We have also ex- 
pressed c(r, r') in terms of its real and imaginary parts as 
c(r, r') = a{r, r') + ib{r, r') and 



Q{r, r', Lu) = du 



(8) 



Equation ([8j is formally exact. The continuum of eigen- 
modes reflects the losses of the PHC cavity. The resonant 
cavity modes arise as sharp resonances in the energy-dependent 
density of these eigenmodes. A very effective approxima- 
tion, close to the resonant mode we are interested in, con- 
sists in writing explicitly its contribution to (O and ac 



LUq 



LUq 



htMC^ 



• Wo 



IQD = 70 + 7r 



(13) 



9" = 



titM 



with 



7r 



the cavity induced damping 



counting for its finite hfetime through a damping constant rate. The quantity c(r, r') models the couphng of the QD 
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Figure 1 Response of the system to a dipolar excitation, 
as given by Eq. |5] (logarithmic scale). In this case uic = 
1062 meV. 



to all other field modes which, in addition to a small en- 
ergy shift, are responsible for the intrinsic radiative lifetime 
of the two-level system, denoted by jr in the CQED for- 
malism. In the expression for g^, the mode volume enters 
through the localized cavity mode wave function <Po{ro). 

With the assumption of a dipolar field being emitted 
close to the dot (virtual oscillating dipole method, see Ref. 
0), i.e. 



£b {rQ,uj) = £o.f{r - ro) 



7 



1 



we get the following result: 



C = 



j dr^l{r)f{r -ro) 



and 

5M = 



fl+ — UJa + ij — (^a + il 



uj — 



uj — Q- 



r(7 



7 - 



(14) 
(15) 

(16) 
(17) 



This is exactly the CQED result |l5]|4l. As seen in Eq. 
( fT4l i. the form factor JF is proportional to the overlap be- 
tween the QD volume and the cavity mode wave function 
through C. Figure [T] displays the response function of the 
system for different dot-cavity detuning D = lOc — ojq, 
showing the typical avoided level crossing. As can be seen 
in Eq. (flJT l. our model provides a very good quantitative ac- 
count of the measured Rabi splitting in realistic situations. 



By using a Gauss-shaped cavity eigenmode, we recover the 
result of Ref. 



9 



(18) 



4 Influence of electron-phonon coupling In a semi- 
conductor, the QD electronic states are coupled to other 
degrees of freedom of the surrounding semiconductor In 
the limit of low excitation and low temperature, coupling 
to acoustic phonons dominates. They can be accounted for 
by using a different expression for x {^)- 




Figure 2 Plot of the quantum dot susceptibility in the 
presence of phonons and without phonons (dashed line), 
T=10K. The cavity mode is also visible (dotted line, for 
u>c = 1061 meV). The scale is logarithmic. 

We consider coupling to longitudinal acoustic (LA) phonons 
via the deformation potential mechanism. This process has 
already been characterized both theoretically ||7]|6| and ex- 
perimentally lis]. Within the second Born approximation 
and restricting to only one phonon band, the QD linear sus- 
ceptibility is given by 161 : 



Xqd (^) 



1 



UJQ 



(19) 



with 



E 



9^ I (1 + ^(9)) 



9q I (Mm)) 



LU + ij — luq — i-ij{q) UJ + — ujQ + uj{q) 



where is the exciton-phonon coupling matrix ele- 
ment. The absorption of the QD (i.e. the imaginary part of 
the QD susceptibility) is plotted in fig.|2] We can observe 
that phonons contribute with a broad background extend- 
ing over a few meV, as compared to the bare QD. In the ab- 
sence of a resonant nanocavity, these sidebands are a van- 
ishingly small feature, difficult to observe in a PL spectrum 
lfT2l . and requiring non linear coherent spectroscopy to be 
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Figure 3 Persistence of the cavity-like peak in the pres- 
ence of phonons (blue line) for different detunings a) D=-l 
meV, b) D=0 meV, c) D=l meV, d) D=2 meV (T=10K). 
The result without phonons is plotted as a dashed line (Wc 
is varying). 




0) (meV) 

Figure 4 Influence of temperature on the phonon induced 
effect (arbitrary units, linear scale), Wc = 1061 raeV . 



characterized fS). In a cavity instead, the highly selective 
spectral filtering of the resonant mode can make this fea- 
ture dominant. In figure [3j we plot the emission spectrum 
of the system for different QD-cavity detunings D and with 
the presence/absence (full line/dashed line) of phonons. 
For zero detuning (panel (b)) both situation show a Rabi 
splitting: in this case, the zero-phonon part of the QD sus- 
ceptibility dominates, and no relevant difference appears. 
For D ~ ImeV, on the other hand, the cavity-like peak is 
strongly enhanced in presence of phonon sidebands. This 
enhancement disappears at still larger detuning, exceed- 
ing the width of the phonon sidebands (see panel d)). This 
proves that phonons can greatly enhance the emission of 
the QD at nonzero detuning. One of the main parameters 
entering this phenomenon is the temperature T. In figure|4] 
we have investigated its influence for a given D ^ 1 meV 
detuning. As expected, the phonon-induced features grow 
with temperature. 

5 Conclusion Many interesting aspects of CQED ap- 
plied to semiconductor systems have been investigated in 
the last decade. However, the scenario in which a QD is 
treated as a macroatom is often oversimplified. In this work, 
we have given a simple example of deviations from the 
usual CQED results, when the microscopic nature of the 
electronic states of the QD is taken into account. We have 
shown that the influence of phonons strongly affects the 
emission properties of a QD embedded in a PHC nanocav- 
ity. Other semiconductor-related mechanisms, like for ex- 
ample Coulomb interaction with the continuum of elec- 
tronic states present in the wetting layer, would lead to 
qualitatively similar effects. All these effects should be thor- 
oughly considered in any study of a QD-cavity system aimed 
at an implementation of quantum information science. 
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